Phenomenological model for ordered onions under shear flow 
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We propose a phenomenological model for the multi-lamellar vesicles (onions) formation induced 
by shear flow. In a nonionic surfactant (C12E4) system, onion phases under a fixed shear flow 
within a certain range show the order-disorder transition accompanied with a size jump by changing 
temperature. Our model can simulate ordered and disordered onion phases with different onion 
sizes. We show numerical results of the onion formation simulated by the model and also discuss 
what factors in this system are critical to cause the transition between these two different onion 
phases. 
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I. INTRODUCTION 

Surfactants in water form various kinds of assemblies 
such as lamellar, micelles and vesicles. Multi-lamellar 
vesicles, which are also called onions, are one of the in- 
teresting morphologies induced by shear flow. The onion 
formation under shear flow has been, in fact, observed 
in experiments [H41]- Onions fill space with various sizes 
and forms of polyhedra. Each onion consists of concen- 
tric lamellar membranes, which were observed by electro- 
microscopy [5]. The size of an onion depends on the 
shear rate and the characteristics of a surfactant mem- 
brane. Several theories have been proposed to estimate 
onion size. Diat et al. proposed the model in which ap- 
plied shear stress is balanced by the elastic stress given 
by the curvature energy of membrane [1] . On the other 
hand, van der Linden et al. proposed the idea that ap- 
plied shear stress is balanced by the energy cost to de- 
form an onion [f|. These models have been examined 
by experiments especially concerning shear-rate depen- 
dence. The formation of onions, i.e. closed-packed or 
disordered structure, is another problem to be investi- 
gated. Panizza et al. suggested that the disordered onion 
phase is made of randomly oriented monodomains which 
consist of closed-packed onions Q ■ 

The characteristics of a surfactant membrane, e.g. 
elasticity, depend on temperature as well as substances. 
An interesting behavior related to the temperature de- 
pendence has been observed in the shear-induced order- 
ing of onions in a nonionic surfactant (C12E4) system: 
Onions form a two-dimensional (2D) close-packed honey- 
comb structure in a certain temperature interval and they 
are disordered below and above the temperature interval 
[H, [^]. Moreover, in recent experiments, it has been re- 
vealed that the order-disorder transition is accompanied 
by a size jump 0. The size of onions is 5-6 times larger 
in the ordered phase than in the disordered one. To our 
knowledge, no theoretical model about the transition has 
been proposed, though the order-disorder transition with 
a size jump had already been reported in other systems 

una. 

In this paper, we propose a model to simulate ordered 
and disordered onion phases with different onion sizes 



from the view point of pattern formation. The model 
provides the means to demonstrate the formation of an 
onion phase in real space. Experiments by the small- 
angle light scattering (SALS) and small-angle X-ray scat- 
tering (SAXS) give only images in reciprocal space. It 
is difficult in those scattering experiments to see what 
structure onions form in three dimensions. Simulations 
in real space will give a clue to analyze the structure of 
the onion phase which is observed experimentally. In this 
paper, we focus on the onion formation induced by shear 
flow. Our model illustrates the onion pattern which is 
stable under shear flow, although the model is not suit- 
able to discuss the lamellar-onion transition. 

From the view point of pattern formation, 2D hexag- 
onal structure such as the ordered onion phase ap- 
pears ubiquitously. One of the most famous examples 
is nonequilibrium fluid dynamics known as Rayleigh- 
Benard convection, of which pattern formation is well 
described by the Swift-Hohenberg (SH) equation. The 
SH equation is a well-known model equation for nonequi- 
librium pattern formation phenomena [ill ], and it has 
a property of so-called potential dynamics. We derive 
a phenomenological model, which is similar to the SH 
equation, of the onion formation by constructing the free 
energy of onion phases under shear flow. 

The rest of this paper is organized as follows. In Sec. 
HI1 we introduce a phenomenological model describing the 
onion formation induced by shear flow. In Sec. IIII1 we 
show numerical results of ordered and disordered onion 
phases with different onion sizes. Section lTVl is devoted to 
the discussion about how the onion size and the ordering 
structure of onions are selected. Finally, we conclude this 
paper in Sec. [V] 



II. MODEL 

We here construct a phenomenological model for the 
shear-induced onion phase. We mainly focus on the tran- 
sition of the onion formation (size and ordering structure) 
and do not consider the transition between lamellar and 
onion phases in our model. To simulate the onion forma- 
tion under shear flow, we derive the free energy of onion 
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FIG. 1: Schematic representation of the shear-induced onion 
phase in our model. 



Free energy for the onion phase under shear 
flow 



On the above assumptions, we give the following Gibbs 
free energy G including the contribution from shear flow. 



G — F a - 



rfacc 



aA. 



(1) 



The free energy G, which is a function of z(x,y), con- 
sists of two parts: (i) the contribution from the free en- 
ergy of a deformed onion (F on i on ) and (ii) the contribu- 
tion from the stress applied to the surface membrane of 
packed onions under shear flow (F sur f aco — a A). We derive 
these two parts of the free energy in what follows. 
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FIG. 2: (Color online) Schematic representation of a distorted 
vesicle surface. Solid and dot-dash curves denote the surface 
of a distorted vesicle and that of a spherical one, respectively. 



phases including the contribution from shear flow, so that 
we can realize the stationary distribution of onions by 
minimizing the free energy. 

We first assume that onions, which have a certain iden- 
tical size R, are packed to form a 2D board. The centers 
of onions are placed on the x-y plane. The spatial dis- 
tribution of onions is expressed as the distribution of the 
height of onion surface, z(x,y) (see Fig. [1]). 

Second, considering onions are closed-packed to form 
a 2D honeycomb structure (Fig. [IJ, we assume that the 
shape of an onion is a slightly-distorted sphere (see Fig. 
[2]) because spherical onions cannot fill space without de- 
formations. In fact, actual onions show a polyhedral 
structure in the experimental observation [5(, and the 
distortion of the membranes is localized at corners. How- 
ever, we assume the distortion of membranes is uniformly 
distributed over slightly-distorted spherical vesicles for 
simplicity. The validity of this assumption is discussed 
in the Sec. [GO 



1. Contribution from a deformed onion 

The free energy of a deformed onion -Fonion consists of 
bending energy fbond and interaction energy F mt coming 
from compression of the layers. We here derive these two 
kinds of free energy as follows. 

We first derive -Fbond in a similar manner to that of Ref. 
The bending free energy per area F/A for a bilayer 
membrane is represented as the following function of the 
mean and Gaussian curvatures, H and H, respectively. 

- = 2kH 2 + RH + -dH* + -c 2 H 2 + 2c 3 H 2 H, (2) 

where k and R are the conventional bending and Gaussian 
curvature moduli, respectively, c\, c 2 and C3 are nonlinear 
moduli. We here consider a vesicle which has a spherical 
shape, which indicates that H 2 and H are equivalent. 
Thus, the free energy, F s h e \\(r), of a spherical shell with 
radius r is given by 

F shcU (r) = 4nkH V + ndH\ 2 , (3) 

where k = 2k + R and c = c\ + c 2 + 8C3. 

We next consider distortion of the spherical shell due 
to packed-onion formation. We assume that the radius 
of curvature changes from an ideal radius r to that for 
a distorted spherical shell r + 5r (Fig. [5]). Since the 
degree of distortion is quite smaller than the ideal radius 
(r 3> Sr), we expand the mean curvature H in a = Sr/r 
and keep terms up to the third order. 

H = — !— ~ -{1-a + a 2 -a 3 ). (4) 
r + Sr r y ' w 

Similarly, we have H 2 and H A as follows. 

H 2 = -^(1 - 2a + ia 2 -4a 3 ), (5) 

H* = 4(1 -4a + 10a 2 - 20a 3 ). (6) 
r 4 

Substituting these expressions into Eq. ([3]), we have 
F shc n(r) = 4ttk(1 - 2a + 3a 2 - 4a 3 ) 

+ — (l-4a+10a 2 -20a 3 ). (7) 
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Thus, the free energy / on ion(-R) (per unit volume of an 
onion with radius R made of n vesicles with layer spacing 
d) is given by 



-FshellfV = jd) 



/onion (R = nd) = ^ 
3=1 

= 37^(1 - 2a + 3a 2 -4a 3 ) 
dR- 

+ 4i&( 1 -^) d-4a + 10a 2 -20a3 

(8) 

where a = SR/ R is the ratio of distortion on the outer- 
most shell of an onion. Let us consider the distortion of 
the onion surface. As illustrated in Fig. [3J the difference 
between the height of the distorted outermost shell and 
that of the spherical one is taken to be e, which leads the 
following relation: R 2 — (z— e) 2 = (R + 5R) 2 — z 2 . Under 
the conditions of slight distortion (R ^> SR and z ^> e), 
we obtain the following relation 



a = SR/R = {e/R 2 )z 



(9) 



By using the above relation, we finally have -Fbend as a 
function of z(x,y), 



Fbc 



dxdy 

= I dxdy 
3c 



/onion (-^) ' Z 



3k 

dR 2 



2az 2 



3a 2 z 3 



4a 3 z 4 ) 



Ad 2 R 3 



1 - — J (z - 4az 2 + 10a 2 z 3 - 20a 3 z 4 ) 

(10) 



We next derive the interaction energy F; n t related to 
compression. The lammellar interaction between layers 
is given by ^B d 2 , where B and d are the compression 
modulus and layer spacing, respectively. Applying the 
interaction to layers inside a onion with radius R, the 
interaction between vesicle layers depends the height of 
the onion surface z(x, y) because effective layer spacing is 
given by j^d. We thus have F int as a function of z(x,y), 



Fint = / dxdy 



n „ / z , 
2 Bo (* d 



(11) 



where we have used n = R/d. 

From the above derivation, we finally obtain the free 
energy of a deformed onion F on ion 



= F 



bend 



= / dxdy 



Fmt 

^J^j(-(4 + 20A)a 3 z 4 + (3 + 10A)a 2 z 3 

(12) 



B 0^~ 2 



-(2 + 4A-_paz 2 + (l + A)z 



where A = izwi and Bo 
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2. Contribution from the stress applied to the surface 
membrane under shear flow 

We next consider the free energy from the stress ap- 
plied to the surface membrane of packed onions under 
shear flow. Suppose that onion surfaces are covered by 
a single membrane as if the membrane is a part of each 
onion surface, and also that the membrane is subjected 
to an effective stress that comes from surface tension and 
(Shear stress. From these assumptions, the free energy 
contribution from the effective stress is written as 



rfacc 



<jA = J dA{2nH 2 - a), 



(13) 



where a and A denote the effective stress and the total 
area of the membrane, respectively. The first term is the 
bending free energy of the onion surface as a single mem- 
brane and the second term is the effective stress applied 
to the membrane under shear flow. We assume the effec- 
tive stress a includes contributions not only from effec- 
tive surface tension but also from shear stress. Although 
the shear can affect vesicle membranes inside onions, we 
impose the effect from shear stress on only the surface 
membrane for simplicity. 

Using s — Vz, we have the following expressions [l3[ 



dA = vl + s 2 dxdy, 
V ■ s 



H = 



(14) 
(15) 



Substituting these expressions into Eq. fjl3|) and assum- 
ing s 2 < 1, we finally obtain 



Surface -aA= \ dxdy 2k(V 2 z) 2 - - (Vz) 2 - a 



(16) 



B. Distribution of the onion surface z(x,y) 

To simulate the stationary distribution of the onion 
surface z(x,y), we use the following Ginzburg-Landau 
(GL) equation (l4T |. in which time evolution is given by 
the variation of the free energy G, 



dz _ 5G _ - 
at dz 



,(z) + (tV 2 z - 4kV 4 z, (17) 



where 

/onion \Z) 



'6 k 



-(16 + 80A)a 3 z 3 + (9 + 30A)a 2 z 2 

(18) 



4d 2 R 3 

- (4 + 8A - B )az + (1 + A) 



Here, the phenomenological coefficient L has the dimen- 
sion of time over mass and we fix L = 1 in the following. 
By solving the above equation, we finally obtain the sta- 
tionary distribution of the onion surface z[x, y). Since 
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TABLE I: System parameters 



a effective stress Bo compression modulus 

k bending modulus R radius of an onion 

R saddle-splay modulus d layer spacing 
5 nonlinear modulus e shift of z caused 

by distortion 



TABLE II: System parameters chosen for large (case A) and 
small (case B) onion phases. 





a 


K 


c 


d 


e 


case A 


20.0 


— 13.0 


-3.9 


0.014 


1.20 


case B 


200.0 


— 1.0 


-3.0 


0.013 


0.35 



the stationary distribution is realized by minimizing the 
free energy G of onion phases under shear flow, we can 
simulate the onion pattern formation under shear flow by 
solving Eq. (IT71) . 

It should be noted that Eq. (fT7|) has essentially the 
same form as the SH equation. In other words, the equa- 
tion we have derived is in the same class as the SH equa- 
tion, although the mechanism of onion pattern formation 
is different from that of the system in which the SH equa- 
tion was originally derived. 



III. RESULTS 

We here show numerical results obtained from our 
model to demonstrate the change between a small- 
disordered onion phase and a large-ordered one. Table 
U lists the system parameters of the model. The system 
size is 200 x 200 (lattice) with the periodic boundary con- 
dition. To integrate the model equation, we have used 
the Crank-Nicolson scheme. Since we are mainly inter- 
ested in the packed-onion ordering, we focus only on the 
stationary distribution of z(x,y) after a sufficiently long 
time. In our numerical simulations, n is fixed (k = 10.0), 
because it is considered to be almost constant under ex- 
perimental conditions for onion phases [15j. Bq and R 
is set by using following relations: B = a 2 d/ k and 
R = 47t-\/2k/3c7. The validity of these relations is dis- 
cussed in the next section. Two different sets of the other 
parameters (cases A and B) are listed in Table HI 

Figure [3] shows the distributions of z{x,y) and spatial 
power spectra for case A and case B. For case A, the size 
of onions is large (R ~ 7.3) and onions show a honeycomb 
structure (panel (a)). By contrast, for case B, the size of 
onions is quite small (R ~ 2.3), although the distribution 
of onions looks an almost perfect honeycomb structure 
(panel (b)). However, spatial power spectra for these 
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FIG. 3: (Upper panels) Spatial distributions of z(x,y) for (a) 
case A, and (b) case B. The brighter the color, the larger the 
value of z. (Lower panels) Spatial power spectra for (c) case 
A and (d) case B. 



cases are totally different. Panel (c) shows a perfect six- 
fold symmetry for case A, although panel (d) shows an 
isotropic circular spectrum for case B. This isotropic ring 
indicates that the onion phase consists of some small do- 
mains of packed onions (i.e. onion texture, @,[Hj]): Each 
small domain has the perfect six-fold symmetry but has 
different directions. These random directions of small 
domains result in the isotropic ring spectrum. 

Furthermore, the time duration before the distribution 
of z(x, y) falls into a stationary state is totally differ- 
ent for cases A and B. While the ordered onion phase 
gradually appears in case A, some small domains ap- 
pear very quickly in case B. The difference can be ex- 
plained in terms of the linear-growth rate, which is es- 
timated from the linear stability analysis. Since the 
linear-growth rate is a very small positive value in case 
A, the stationary pattern can reach the potential min- 
imum, where the ordered honeycomb structure is sta- 
ble. On the other hand, in case B, because the linear- 
growth rate is a relatively large positive value, the sta- 
tionary pattern quickly reaches a metastable state, where 
some randomly-directed monodomains coexist, and stays 
there. 

We have also checked the robustness of these results by 
simulations in a larger system size. These results indicate 
that our model is able to reproduce the change between 
the large-ordered onion phase and the small-disordered 
one by changing the system parameters. We next discuss 
the validity of changing these parameters. 
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IV. DISCUSSION 

We here discuss what factors in our model are essen- 
tial to induce the change between the large-ordered onion 
phase and the small-disordered one, focusing on tha va- 
lidity and temperature dependence of the parameters in 
our model. 



A. Onion size 

At first, we focus on how to determine onion size R in 
our model. As mentioned above, Eq. (|17|) has a simi- 
lar mathematical structure to that of the SH equation, 
which is a model equation to simulate periodic domain 
patterns. The characteristic length Ao of the domain pat- 
tern produced by the SH equation can be estimated ef- 
fectively from the li near s tability analysis. In the case of 
Eq. ([T7| . Ao = Ao represents the distance be- 

tween the centers of neighboring onions and it also has a 
relation Ao = \J~ZR because of the geometrical structure 
of packed onions (see Fig. [!}. Therefore, the estimated 
onion size is given by 

This equation indicates that the onion size in our model 
is determined by the ratio of k and a. We set R in our 
simulations by using Eq. (TT9l) . 

Actually, only a is the essential parameter to determine 
the onion size, because n is fixed in our simulations. Un- 
der a fixed shear rate, it is assumed that the change of a 
is caused by the change of the surface tension. Accord- 
ing to experimental studies @,[I3|; * ne elastic modulus G' 
measured in linear viscoelastic experiments is associated 
with the surface tension er s and onion size R: G' ~ a s /R. 
Assuming a s ~ a and substituting Eq. (fl"9| into the 

3 

above relation, we have G' ~ er| . In the experiment of 
the C12E4 system, on the other hand, G' decreases as 
the radius of onions increases (l5j | . These experimental 
results indicate that a s in the large-ordered onion phase is 
smaller than that in the small-disordered one. This prop- 
erty of the surface tension <j s corresponds qualitatively to 
that of the parameter a in our simulations (a = 20.0 and 
200.0 in cases A and B, respectively). 

We further discuss the change in onion size as a func- 
tion of shear rate 7 and layer spacing d. Diat et al. es- 
timated onion size R by considering a balance between 
the elastic stress and the viscous stress as follows [ij. 

„ / 4ttk 
V rjdj 

where 7/ is the solvent viscosity. For a fixed 77 (i.e. New- 
tonian fluid), this relation implies R ~ (dj)~z. This 
estimation accords with 7" 2 dependence experimentally 
observed by themselves. Assuming the effective stress a 
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FIG. 4: Phase diagram of onion phase as a function of onion 
size R and layer spacing d. The solid curve indicates the 
boundary between the onion phase ("Onion") and a phase 
where no pattern appears ("No pattern"), which is calculated 
from the linear stability analysis. The slope of the dashed line 
is -§. 



in our model is estimated as a ~ 777 and substituting it 
into Eq. (flU)) . we also have R ~ 7" 2 for a fixed 77. We 
further estimate the relation between onion size R and 
layer spacing d from the linear stability analysis of Eq. 
dTTJ) . Figure 0] shows a phase diagram of onion phase as 
a function of R and d. The boundary of onion phase, 
which is calculated from the linear stability analysis, has 
a slope — \. It indicates the relation R ~ because 
the ordered onion phase appear near the boundary, that 
is, the region where the linear-growth rate is very small. 

Actually, the order-disorder transition with a size jump 
is related to the compression of the onion lay ers. It causes 
a discontinuous change of the viscosity [18| . The experi- 
mental result of Ref. [TH also shows the shear thinning, 
that is, the viscosity 77 depends on the shear rate 7. How- 
ever, the effective stress a in our model includes not only 
the shear stress but also 77, 7, and d. Thus, we cannot 
discuss the compression or shear thinning in our model. 



B. Degree of the distortion 

We next consider the validity of the distortion of vesicle 
membranes due to packing of onions. Although we as- 
sume that membranes of packed onions are slightly dis- 
torted from a sphere, actual onions have a polyhedral 
structure [5J . We clarify the validity of the assumption 
by estimating the degree of the distortion e/R. Inside an 
onion formed by polyhedral vesicles, the effective surface 
tension <r e ff and the penetration length £, which corre- 
sponds to the distortion of vesicle membr anes, can be 
estimated as a c ft — \/KB and £ ~ y/KjB, respectively, 
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where K is the bending modulus K — n/d and B is the 
compression modulus @,LL9|]. We here assume that the ef- 
fective surface tension corresponds approximately to the 
stress a in our model, namely a c g ~ a. We thus have 
£ ~ njad. In addition, we can estimate the compres- 
sion modulus Bq in our simulation by using the above 
estimation of a c g, and consequently, we set Bq = a 2 d/n. 

On the other hand, the onion size R can be estimated 
to be R ~ \/n/a as mentioned above. We thus obtain 
the following relation: ~ R. This indicates that 

the larger the onion size R the larger the distortion e, if 
one assumes £ ~ e. The above relation is also equivalent 
to the relation £/i? 2 = constant. This corresponds to 
our assumption a — e/R 2 = constant. Therefore, the 
assumption for the distortion of vesicle membranes due 
to packing of onions is consistent with other assumptions. 
In addition, we set e in our simulations to satisfy the 
above relation that the larger onion has larger distortion. 



C. Elastic properties of vesicle membranes in 
packed onions 

Elastic properties of vesicle membranes in packed 
onions, such and c, should also be considered. 

We first focus on the signs of coefficients of the bending 
free energy such as n and R. In general, the saddle-splay 
modulus R is a negative value unless the sponge (L3) 
phase is stable and spherical vesicles are stabilized for 
k = 2k + R < 0. However, the effects of shear flow are 
not taken into account in this condition. Thus, it is not 
necessary in our model to satisfy the condition. We have 
actually checked simulations for different sets of param- 
eters k, R, and c and consequently found that the onion 
phase can be obtained only when R < and k > 0. This 
result suggests that the shear stress stabilizes unstable 
spherical vesicles for k > 0. 

We further refer two important nonlinear effects of 
the bending free energy to simulate onion phases in our 
model: nonlinear modulus in the bending energy and 
the distortion of membranes due to packing of onions. 
The usual expression of the Helfrich curvature free energy 
does not include the nonlinear coefficient c. However, c 
should be considered in an onion phase, because high 
curvatures are expected near the core of an onion [12l |. 
Also in our model, c is one of the essential parameters to 
simulate onion phases, because no onion phase is gener- 
ated from our model if c is taken to be 0. With regard to 
the distortion of membranes, it is also necessary to con- 
sider the distortion of vesicle membranes in our model 
to simulate onion phases. The third-order expansion of 
the curvature on a distorted vesicle is necessary to gener- 
ate the ordered onion formation, because the non-linear 
terms of z in Eq. (fT8|) play a similar role of a double- well 
potential. 



D. Temperature dependence of the parameters 

We further consider the relation between the tempera- 
ture dependence of the system parameters in our simula- 
tions and that of the physical quantities in experimental 
studies. Since a is the critical parameter to determine 
the onion size in our simulations, it should have a simi- 
lar temperature dependence to that of the onion size. In 
other words, a should be minimized in a certain temper- 
ature interval and have larger values above and below the 
interval. 

On the other hand, Kosaka et al [13] reported that 
layer spacing d can change around the lamellar-onion 
transition, and also discussed temperature dependence 
of R. They measured d at different concentrations of the 
surfactant and showed that the ratio S/d of d and the 
thickness of bilayers 8 has a temperature dependence. 
Although they focused on the lamellar-onion transition 
in their experiments, the temperature dependence oiS/d 
may be one of the important factors to cause the order- 
disorder transition in onion phases. We therefore change 
the value of R and d in our simulations, assuming the tem- 
perature dependence of these parameters. Actually, they 
are not essential parameters to determine the onion size, 
but important parameters to determine the linear-growth 
rate, that is, determine whether onion phases appear or 
not. 



V. CONCLUSION 

We have proposed the phenomenological model for the 
onion pattern formation induced by shear flow. We have 
derived the free energy of an onion phase by consider- 
ing contributions from deformed packed onions and from 
the stress applied to the surface membrane on packed 
onions under shear flow. Our model includes two im- 
portant nonlinear effects, that is, nonlinear modulus c 
in the bending energy and the distortion of membranes 
due to packing of onions. These nonlinearity are crucial 
in simulations to reproduce onion phases in our model. 
Then, we have demonstrated that the two different onion 
phases, large-ordered one and small-disordered one, can 
be simulated by adjusting the system parameters in our 
model. The ordered onion phase has a clean honeycomb 
structure made of packed onions. The disordered onion 
phase consists of small domains of ordered packed onions. 
Each domain has the six-fold symmetry, and domains are 
laid in random directions. The size of packed onions is 
determined by the ratio of the bending modulus k and 
the effective stress a. Thus, a is the main control param- 
eter under the fixed-K condition, though other parame- 
ters, such as R and d, are also important for emergence 
of onion phases. In fact, the system parameters in our 
model, such as cr, are difficult to compare with the physi- 
cal quantities which can be measured in experiments. In 
that sense, our model can be improved, nevertheless, it 
is useful as a simple model for the shear-induced onion 
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formation. 
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